home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cspr.z / cspr
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCSSSSPPPPRRRR((((3333FFFF))))                                                              CCCCSSSSPPPPRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CSPR - perform the symmetric rank 1 operation   A := alpha*x*conjg( x' )
  10.      + A,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CSPR( UPLO, N, ALPHA, X, INCX, AP )
  14.  
  15.          CHARACTER    UPLO
  16.  
  17.          INTEGER      INCX, N
  18.  
  19.          COMPLEX      ALPHA
  20.  
  21.          COMPLEX      AP( * ), X( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      CSPR    performs the symmetric rank 1 operation
  25.  
  26.      where alpha is a complex scalar, x is an n element vector and A is an n
  27.      by n symmetric matrix, supplied in packed form.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      UPLO   - CHARACTER*1
  32.             On entry, UPLO specifies whether the upper or lower triangular
  33.             part of the matrix A is supplied in the packed array AP as
  34.             follows:
  35.  
  36.             UPLO = 'U' or 'u'   The upper triangular part of A is supplied in
  37.             AP.
  38.  
  39.             UPLO = 'L' or 'l'   The lower triangular part of A is supplied in
  40.             AP.
  41.  
  42.             Unchanged on exit.
  43.  
  44.      N      - INTEGER
  45.             On entry, N specifies the order of the matrix A.  N must be at
  46.             least zero.  Unchanged on exit.
  47.  
  48.      ALPHA  - COMPLEX
  49.             On entry, ALPHA specifies the scalar alpha.  Unchanged on exit.
  50.  
  51.      X      - COMPLEX array, dimension at least
  52.             ( 1 + ( N - 1 )*abs( INCX ) ).  Before entry, the incremented
  53.             array X must contain the N- element vector x.  Unchanged on exit.
  54.  
  55.      INCX   - INTEGER
  56.             On entry, INCX specifies the increment for the elements of X. INCX
  57.             must not be zero.  Unchanged on exit.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCSSSSPPPPRRRR((((3333FFFF))))                                                              CCCCSSSSPPPPRRRR((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      AP     - COMPLEX array, dimension at least
  75.             ( ( N*( N + 1 ) )/2 ).  Before entry, with  UPLO = 'U' or 'u', the
  76.             array AP must contain the upper triangular part of the symmetric
  77.             matrix packed sequentially, column by column, so that AP( 1 )
  78.             contains a( 1, 1 ), AP( 2 ) and AP( 3 ) contain a( 1, 2 ) and a(
  79.             2, 2 ) respectively, and so on. On exit, the array AP is
  80.             overwritten by the upper triangular part of the updated matrix.
  81.             Before entry, with UPLO = 'L' or 'l', the array AP must contain
  82.             the lower triangular part of the symmetric matrix packed
  83.             sequentially, column by column, so that AP( 1 ) contains a( 1, 1
  84.             ), AP( 2 ) and AP( 3 ) contain a( 2, 1 ) and a( 3, 1 )
  85.             respectively, and so on. On exit, the array AP is overwritten by
  86.             the lower triangular part of the updated matrix.  Note that the
  87.             imaginary parts of the diagonal elements need not be set, they are
  88.             assumed to be zero, and on exit they are set to zero.
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.